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Abstract: To properly discover a charged Higgs Boson (H^) requires its spin and cou- 
plings to be determined. We investigate how to utilize ti spin correlations to analyze the 
couplings in the decay t — > bH + — * br + v T . Within the framework of a general Two- 
Higgs-Doublet Model, we obtain results on the spin analyzing coefficients for this decay 
and study in detail its spin phenomenology, focusing on the limits of large and small values 
for tan (3. Using a Monte Carlo approach to simulate full hadron-level events, we evaluate 
systematically how the — * r^v T decay mode can be used for spin analysis. The most 
promising observables are obtained from azimuthal angle correlations in the transverse rest 
frames of t(t). This method is particularly useful for determining the coupling structure of 
in the large tan/3 limit, where differences from the SM are most significant. 
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1. Introduction 

Finding a fundamental spin zero boson with electric charge would be a direct sign of physics 
beyond the Standard Model (SM). The existence of a charged Higgs boson pair (i/ ± ) with 
these properties is predicted by Two-Higgs-Doublet Model (2HDM) extensions of the SM 
Higgs sector. The primary motivation for studying the 2HDM is supersymmetry, which 
requires an even number of Higgs doublets for cancellation of triangle anomalies. Charged 
Higgs searches at hadron colliders are divided into two regimes, separated by the dominant 
mode of production. When H is heavy (mjj+ > m-t), it is produced primarily through 
the gg — > H + ib and gb — > H + t processes [1-4]. When, on the other hand, H is light 
(m H + < mt — mb) and the decay t — > bH + [5-9] opens up, this quickly becomes the 
dominant production mode. 

The most stringent model-independent limit on the mass of from a direct search 
experiment comes from LEP: m^± > 79.3 GeV [10] at 95% CL, assuming only the decays 
H + — > cs and H + — > t + v t are possible. Even tighter constraints on m H + have later 
been derived using Tevatron data [11], but these are not independent of the other 2HDM 
parameters. Neither are indirect constraints on obtained from i?-physics observables. 
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The upcoming searches for H planned by the LHC experiments will have good sen- 
sitivity to discover over a wide parameter range [12-16], especially when is light. 
However, even if some candidate state was to be found, this discovery alone would not 
be enough to establish the validity of the Higgs mechanism as described by the 2HDM. To 
do this requires further that the spin, and the couplings, of this new particle be determined. 
Here we investigate one possibility to address this issue when m#+ < mt — m-b- In addition 
to providing a discovery channel, the t — ► bH + decay mode could modify the ordinary V-A 
Lorentz structure of weak top decay significantly. As we will show, this fact can provide a 
handle on the spin and coupling structure of H by making use of spin correlations. 

Top quarks produced in pairs at hadron colliders constitute an interesting laboratory 
for observing spin effects in high-energy physics. Since the timescale for weak top decay 
1/Tt is much shorter than the typical hadronization timescale 1/Aqcdj the heavy quarks 
will decay before hadrons can form [17]. No hadronic effects will therefore obfuscate the spin 
information. Unlike the case for the lighter quarks, this fact allows for reliable perturbative 
calculations of the relevant spin observables. Furthermore, since the charged current weak 
interaction violates parity maximally, the decay self-analyzes the spin of the top quark. 
This means the full spin information will be imprinted in the angular distributions of the 
different decay products. 

In order for the angular information to be useful for investigating the couplings in- 
volved, a method to determine the spin projection of the decaying top quark is required. 
At a hadron collider, this can be achieved by exploiting correlations between the top quark 
spins. By using the decay information from one side of a ti event, it is possible to deter- 
mine statistically the polarization of the other top. Spin correlations in top quark pair 
production and decay have been extensively discussed within the Standard Model (SM) for 
hadron-hadron colliders [18-30], and for e + e~ experiments [31-36]. Utilizing the top spin 
information to study physics beyond the Standard Model was considered in the context 
of anomalous VTt6-couplings [31,37,38], for probing extended Higgs sectors and effects of 
CP-violation [39-44] , and within theories with extra dimensions [45,46]. 

We study the spin phenomenology of a light scalar sector in full generality, ignoring 
indirect constraints on the 2HDM. We will however restrict the mass to m#+ > mw, 
as required by the non-observation of H in direct search experiments. Due to the, in 
the context of spin observables, relatively limited sample of ti events available from the 
Tevatron, our main focus will be on prospects for observations at the LHC. We try to 
comment on issues that are of relevance also for the analysis of Tevatron data. 

The organization of this paper is as follows. First, in Section |2[ we discuss ti pair 
production at hadron colliders, and how spin correlations come about in this process. 
Then in Section |3] we introduce the phenomenological model, followed by a brief review 
of the relevant theory for polarized top quark decay, and results on the spin analyzing 
efficiencies in models where the top quark can decay through a charged Higgs boson. 
Section |I| describes a Monte Carlo simulation study of these effects, and discusses possible 
experimental observables. Finally, Section |5| contains a summary and the conclusions of 
this work. 
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2. Top Quark Pair Production at Hadron Colliders 



We will adopt the latest combined value for the top mass = 170.9 ± 1.8 GeV [47]. Pair 
production of top quarks occurs in leading order QCD both through qq annihilation via 
an s-channel gluon, and through gluon fusion for which s-, t- and u-channel exchanges are 
possible. The total hadronic cross section has long been known to NLO accuracy [48,49]. 
With CTEQ6 parton distributions [50] and a common choice of scales /ur = = frit, 
we obtain a(pp — > ti + X) ~ 900 pb at yfs = 14 TeV using the NLO MC generator 
POWHEG [51]. One year of LHC running at low luminosity, corresponding to J C = 10 
fb _1 , will therefore produce of order 10 7 ti events. 

Individually unpolarized top quarks are still produced with strong correlations between 
their spin projections in a suitable basis. The nature and magnitude of these correlations 
depend on the partonic center of mass (CM) energy. Following [24] we define the production 
correlation, as a function of the invariant mass M t i of the top pair, to be 



Cij(M?r) = ^¥j + tih)-^(hh + ¥ } ) (21) 



for the partonic subprocess involving initial state partons Arrows indicate the spin 

projection on the chosen spin quantization axes (which may be different for t and t). We 
choose here to work exclusively in the helicity basis, in which the spin is quantized along 
the momentum directions of the t(t) in the partonic CM frame. In this basis, the notation 
(R,L) is sometimes used interchangeably with (j, j) to denote the two spin projections. 

Near threshold, the ti are always produced in an S-wave state. For production dom- 
inated by qq — > ti through an s-channel gluon, the overall angular momentum state will 
therefore be 3 Si. Out of the three states composing the triplet, two correspond to opposite 
helicities for the two top quarks, whereas one state gives equal helicities. The combined 
correlation according to (2J.) is thus C q q{^m^) = —1/3 at threshold. When instead gg — > ti 
production dominates, which is the case for LHC energies, the situation at threshold is re- 
versed. The top pair is now produced in a singlet 1 So configuration, since the initial state 
gluons do not populate J = 1 states [20]. In this case the ti must always come with the 
same helicities, which means that C gg {Am1) = 1. Finally, in the ultra-relativistic limit 
(M^t S> 4mf), helicity conservation requires the ti to have opposite helicities independent 
of partonic subprocess. Hence the correlation Cij(M?r) — ► — 1 for all in the high 

energy limit. 

In general, the total statistical correlation C in a sample of ti events is obtained from 
averaging over the invariant mass of the top pair, while parton distribution functions de- 
termine the relative contributions of the competing production processes: 

C(s) = — Y] / dx 1 dx 2 [a ij (t^ +t l i l ) - aijit^ + t^i l )]f i (x 1 ,fil)f j (x 2 ,fJ>F). 
i,3={q,q,9} 

(2.2) 

The Cij for all partonic subprocesses have previously been calculated to NLO in QCD. 
From these one obtains C = 0.326 [29] in the helicity basis for pp collisions at y/s = 14 
TeV. The residual uncertainty in this number from PDF and scale choices is of order one 
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Figure 1: The dashed red curve shows the resulting value of C in the helicity basis, at the LHC, 
when a cut M t j < M cut has been applied. The black curve shows the fraction of the total cross 
section which survives the cut. 

percent, a value similar to the difference from the LO calculation which gives C = 0.319. 
For Tevatron run-II conditions (pp collisions at y/s = 1.96 TeV) the same helicity basis 
correlation becomes C = —0.352. At Tevatron energies, where qq annihilation totally 
dominates ti production, there exist also more efficient bases for spin quantization in which 
correlations as large as C = 0.8 can be obtained [27]. 

As suggested by the discussion above, it can be beneficial to introduce an experimental 
cut on M t i to increase the spin-purity of the ti sample at the cost of decreased efficiency [23] . 
To increase the component of like-sign helicities in the gluon sample at the LHC requires a 
cut on the maximum M t j. The effect of such a cut on the correlation parameter C is shown 
in Figure [l]. To illustrate the trade-off between spin-purity and efficiency, Figure Q also 
presents the fraction of the total cross section which passes a cut on the maximum M t i. 
As we have already indicated, ti statistics will not be the limiting factor at the LHC. It is 
therefore good to keep in mind that C can be increased using this technique, although we 
do not make explicit use of this fact here. 

3. Spin Information in Top Quark Decay 

In the SM, the top quark decays almost exclusively via the charged current V-A vertex 

1 5 

C wtb = ^V tb W+i^^-b + h.c, (3.1) 

where Vtb ~ 1 is the appropriate element of the CKM matrix. In addition to the Standard 
Model decay, the possibility exists that the top quark decays anomalously. The decay 
could then contain a small V+A component, or it could be mediated by additional bosons 
of different spin. The study of spin correlations opens a window on both these possibilities. 
Our aim here is to explore the latter case, allowing for an extended scalar sector. 

3.1 Charged Higgs Model 

Introducing a charged scalar pair H^, their interactions with fermions are parametrized 



- 4 - 



Table 1: CP-invariant couplings of charged Higgs bosons to fermions in 2HDM (I), 2HDM (II) 
and SUSY-improved type II model (II). 



Coupling 


2HDM (I) 


2HDM (II) 


2HDM (II) 


A 


m u cot (5 


m u cot P 


m u cot P [1 — e' t tan P] 


B 


—rrid cot P 


nid tan P 


m d tan fi 
l+tb tan fi 


C 


mi cot P 


mi tan p 


mi tan P 



by an effective Lagrangian density 



C H = - 9 J V [A (1 - 75 ) + B(l +T B )]d + H-d\B*(l- 75) + ^*(l + 7 5 ) 

2V2m ^ {u,c,t} { 1 J L 

{<M,b} 

[™i +75)i+ *-c-i(i- 7 «,4 

(3.2) 

The A, B, C, and their complex conjugates, are in principle free parameters determining 
the Lorentz structure of the couplings. Note that this model does not assign definite parity 
to the unless B = ±A. For AB = parity is violated maximally. Assuming CP- 
invariance of the scalar sector, the coupling parameters can all be taken as real numbers. 

A model such as ([T^) occurs, for example, as the charged Higgs-fermion sector of a 
two Higgs Doublet Model (2HDM). Here, the SM Higgs sector is augmented with another 
complex Higgs doublet, resulting in two charge conjugate (H^) and three neutral (h, H, A) 
states occurring as physical bosons. To ensure the augmented SM does not allow for tree- 
level FCNC's, certain restrictions apply on how to couple the extended Higgs sector to the 
fermions. For our purposes, it suffices to say that two options are generally considered: 
In the so-called type I model [2HDM (I)], only one doublet is coupled directly to the 
fermions. In the type II model [2HDM (II)], one doublet is coupled only to up-type fermions, 
whereas the other doublet couples only to down-type fermions. The number of independent 
parameters in the Higgs sector is thereby restricted to two at leading order. We will adopt 
for these the ratio tan/3 = V2/V1 of the two doublets vacuum expectation values, and the 
charged Higgs mass mjj+. For the two model types, the charged Higgs-fermion couplings 
are then given in Table |]. Mass parameters appearing in the couplings should be evaluated 
at a scale Q 2 = m 2 H+ , using the MS masses m g (Q 2 ) to ensure proper resummation of large 
logarithmic QCD vertex corrections [52-54]. Since the couples proportionally to the 
fermion mass, we will only be concerned with third generation particles in the following. 

As discussed in the introduction, one possible extension of the Standard Model where 
a 2HDM occurs naturally is the Minimal Super symmetric Standard Model (MSSM). The 
MSSM contains a 2HDM (II), but since the supersymmetry introduces additional particles, 
the two-parameter picture of the 2HDM (II) works only as an effective tree-level description. 
It has been shown [52,55], that quantum corrections due to SUSY-QCD loops can be quite 
sizable for large values of tan p. This holds even in the decoupling limit when all SUSY 
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masses are taken to infinity. How these corrections enter into the effective couplings 
can be seen from the third column of Table [l|. We will call the 2HDM which includes the 
tan/3 enhanced SUSY corrections the "modified type II", or simply II. The corrections are 
of two types: first the so-called eb correction to the relation between the bottom quark mass 
nib and the bottom Yukawa coupling yb- It is caused by gluino-sbottom and chargino-stop 
loops. At one-loop, the dominant contributions to this correction are given by [55, 56] 



2a s n 
3-7T m„ 



H 2 



nir nir 

01 02 



nig m g 



Vt 
16vr 2 



U a2 —H 2 



A a 



ni-+ 

> Xa 



rn 



Xa , 



V a2 , 



(3.3) 



which introduces a dependence on the trilinear coupling A t , the top Yukawa coupling y t , 
and the \i parameter from the superpotential - in addition to the dependence on several 
of the sparticle masses. The real matrices U and V diagonalize the chargino mass matrix. 
The function H 2 is given by 

x\nx 



H 2 (x,y) 



+ 



ylny 



(l-x){x-y) (l-y){y-x)' 



(3.4) 



In the limit when all SUSY parameters and sparticle masses are of similar scale Msusy 
one obtains |e&| ~ a s (Q = Msusy) / (Sir) ~ 10 -2 . The sign of e& is determined by the sign 
of [1. 

The second contribution which modifies the couplings is [55] 



2a s \x 
3ir m g 



t b 



+ 44 H * 



vl m A b 



IQtt 2 4a 



An 



t b 




+ 44 H 2 




+ &! H 2 



Na3, 



(3.5) 



sin 9g and Cq 



where the matrix iV diagonalizes the neutralino mass matrix, Sq 
for the squark mixing angles 9g. The squark mass eigenstates are given by q± = Cqq^ + Sq-q^ 
and q 2 = —Sq-qi, + c^r, with > niq 2 . We note that e' t is numerically similar to e& in 
the case with a common scale for the SUSY parameters. In Section |3.5| , we will return to 
these SUSY corrections when discussing numerical results for the 2HDM. As we will show, 
it turns out that their effects on ti spin correlation observables are small. 



3.2 Top Quark Decay with Polarization 

Assuming that the full width Yt of the top quark, including the scalar decay mode, is 
still very small (Tt/ m t ^ 0.01) we use the narrow width approximation to factorize the 
production from the decay of the heavy quarks. The branching fractions for t — > bH + 
are shown in Figure ^ for 2HDM type (I) and (II). It is clear that type (II) is interesting 
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Figure 2: Branching fractions for t -> bH + in 2HDM (I) (left) and 2HDM (II) (right), as a function 
of tan/3. The contours show m H + = 80 GeV (short-dashed), m#+ = 120 GeV (long-dashed), and 
m H + = 160 GeV (solid). 
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Figure 3: Notation for momenta in semi-leptonic top quark decay mediated either by W + or 

In hadronic decay of the boson, we use the same notation while making the replacements of weak 

isospin partners: l + — > J, and i/j — * u. 



both for small and large tan j3 values, whereas the 2HDM (I) only allows a significant 
B1Z(t — > bH + ) for small tan/3. In the following we will mostly be concerned with the type 
(II) model. 

We treat the decaying ft as independent decays through well-defined channels without 
interference effects. Strictly speaking, a more complete formalism involving off-diagonal 
W ± /H ± propagator elements could be used when mjj+ ~ mw- Ignoring such complica- 
tions, the full structure of density matrices in the 2 — ► 6 matrix element becomes 

\M(2 - 6)| 2 = [R X y^(2 -> ft)® ^(t 3)®^,(t^3)]. (3.6) 



i? is here the fully helicity-dependent spin density matrix for it production. The p(p) are 
decay density matrices of t(t), where i,j label the available decay channels and A, A' (k, k') 
are helicity indices for the t(t). 

To obtain the decay density matrix for a given channel, we use the techniques described 
in Appendix [A]. With momenta defined in Figure [||, the leading order decay density matrix 
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elements for semi-leptonic weak decay of the top quark are given by 



U w ,, , TT -+ ,,+ .,2 2 9w\ V tb\ 2 (P ■ k 2)(h ■ h) . , ,„ „ 

P\\> = \M X y(t^bW + ^bl+v^l = " S xy + k 2 a xx , (3.7) 

W nb W> ' w w 



when the spins of all outgoing particles are summed over. The unit 3- vector k 2 is given 
in the rest frame of the decaying quark. For hadronic decay of the W boson, the matrix 
elements are exactly the same if a) all final state masses are neglected and b) the leptons 
are replaced by their quark counterparts in terms of weak isospin. CP-invariance of the 
decay ensures that \M\\>(t — * bl + V[)\ 2 = \My\(i — * bl~v{)\ 2 . 

Reckon similarly the elements of the decay density matrix when the decay is mediated 
by a charged scalar as defined by the model Q3.2|). In this case, the elements become 



AB 45 
+ A 2 + B 2 1 - £ 



(q 2 - m 2 H+ ) 2 + m 2 H+ T 2 H+ 2m^ 

x A 2 -B 2 / AB 45 VSa a 1 (:LS) 
5 ™-JP + &{ 1+ A* + B»l-t) kiaxx ' 



Here the notation £ = m 2 H+ /m 2 and 5 = m^/mt is used. Let us also introduce a convenient 
short-hand 

for the threshold factor. This function has the general properties /(£, A,B) = 1 for m H + -C 
m t , and /(£, A, B) -> for £^1, unless AB = 0. 

3.3 W Boson helicity 

The perhaps most direct test of V-A theory in top quark decay is offered by examining 
the polarization states of the W boson mediating the decay. Due to the large Yukawa 
coupling yt ~ 1, a fraction m 2 /{m 2 + 2m 2 v ) ~ 0.69 of the W bosons are expected to be 
longitudinally polarized, while the remainder carries a left-handed helicity in the t rest 
frame. Uncertainties in these numbers from higher order corrections, including virtual 
2HDM and SUSY effects, are under control at the 1% level [57,58]. 

When the W decays further, the angular dependencies of the decay products on the 
different helicity states are given by the Wigner (/-functions for the spin 1 representation. 
Combining this knowledge with the polarized matrix element for t — > bW + , the normalized 
lepton angular distribution in leptonic decay of the W is given by 

1 dN(W ^ lui) 3 , 2 2 2 21 / Qin \ 

= — — ^ ttt \ m t sin 0, + Mw(l — cos 0, I , (3.10) 

N dcosflf 4{m 2 t +2M 2 v ) V 1 1 WK /;J ' V ' 

where 0^ is defined in the W rest system as the angle of the lepton momentum to the W 
helicity axis. Using the fact that, in the rest frame of the decaying top, the recoiling b 
quark has its momentum anti-parallel to that of the W, the lepton helicity angle cos 6^ can 
be determined by the invariant product [59] 

cosd l=T — 7 1 ~i~ 7 \ 3.11 

h ■ (fe 2 + fe 3 ) 
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Figure 4: Angular distribution of r leptons in the rest system of the boson mediating top quark 
decay. The curves show the SM (black solid), and SM+2HDM (II) with m H + = 80 GeV for three 
different values of tan/3: 50 (dotted), 25 (short-dashed), and 10 (long-dashed). B1Z(H + — > t + v t ) = 
1 was assumed for all values of tan [3. The contours arc equivalent in the low tan (3 regime, whenever 
BTZ{t ->■ bH+) x BTZ(H+ ->t+v t ) is the same. 



if the b mass is neglected. Assuming further that the decay is mediated through an on-shell 
W, the approximate expression 

cos 9? ~ — 5 5 1 (3-12) 

™>i - m w 



can be obtained from the kinematics of the decay. The form ( 3.12| ) is experimentally 
advantageous since no knowledge of the neutrino momenta is required to determine cos 9[ . 
Being only an approximate on-shell relation, the values obtained using this expression may 
in reality be such that | cos# ; *| > 1 for some events. 

In the decay of a charged Higgs boson, the decay products should be isotropically 
distributed in cos#*. This offers a clear signature for a new charged boson to have spin 0. 
However, even with a large branching ratio t — > bH + , this appreciable difference would 
not contribute much to measurements of the distribution (|3. 10 ) using electrons or muons, 



simply because in the tan/3 regions of interest, BTZ^ — > t ~ 1 in the 2HDM. 1 
If, at the LHC, evidence starts to gather in favor of a light if ± , it would therefore be 
interesting to study the angular distribution of r leptons exclusively using the hadronic r 
decay. Experimentally this presents a formidable task, since the presence of two final state 
neutrinos in the r channel introduces ambiguities in the reconstruction of the r momentum. 
Furthermore, when m#+ > mw, the kinematic assumptions behind Equation (3.12) are 



no longer valid. It is then natural to exploit these kinematic differences fully and treat the 
two decays separately also in the angular analysis in order to establish the spin nature 
of the presumptive H^. 

In Figure we show the angular distribution of r leptons in the rest system of the 
boson mediating the t(t) decay. We show here the expectations for the SM, given by 



Equation 3.1C , and mixtures of SM+2HDM (II) with B1Z(t — ► bH + ) corresponding to 



1 A small contamination from H — > r v T — > l^viv T v T would of course be present also in the lepton 
samples in these cases. 
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three different values of tan (3. It is assumed either that m#+ = 80 GeV, or that kinematic 
effects can be compensated for on event-by-event basis. Even for B1Z(t — > bH + ) < 0.1, 
the events give a significant contribution since BTZ{W ± — ► r ± ^ T ) = 0.1125 is equally 
small. 

3.4 Polarization Observables 



From the matrix elements ( |3.7D and ([3^) some useful hints are obtained on how the spin 
is analyzed in top quark decays. The spin-dependent term, proportional to aw, appears 
with different associated momentum directions in the two channels. The particle with this 
momentum will analyze the spin most effectively in the corresponding case. To see how 
this comes about, recall [60] how the spin of a given top quark is analyzed. The decay 
products will experience angular distributions reflecting the spin state of the parent. Each 
polarized partial width T(t^ — > bX + — ► bl + vi) of a decaying fermion can be put in the form 

1 dr l + aiCos6*j 

rw; = 2 ' ( } 

where 9{ denotes the angle of decay product i to the spin helicity axis, calculated in the rest 
frame of the decaying particle. The spin analyzing coefficients on determine the efficiency 
of a given particle to analyze the spin of the parent. The factorization of T in energy- 
dependent and angular parts holds to a high degree also when including radiative QCD 
corrections [60]. 

To obtain the full set of on for a given decay, it is necessary to integrate the polarized 
matrix elements. Using the kinematic variables x = 2p ■ kijvn 2 and y = (&2 + k^) 2 /m 2 , the 
Dalitz parametrization [61] of the 3-body phase space is written as 

d$ 3 = — — dxdyd7d/3dcos6». (3.14) 
(27r) t) 

The three Euler angles 7, 0, and 9 are here chosen according to [62], so that 9 coincides 
with the helicity angle 9i discussed above. The integration over is always trivial and 
gives 2ir. The integration over 7 is non-trivial only for two, spin-dependent, quantities. 
With positive spin projection along the helicity axis, the result is 



d7 s ■ ki = 2tt cos 9 



m. 



and since Yli^i ■ s = p ■ s = 0, with k-z ■ s = — cos9^m 2 , the other interesting integral 
becomes 



/ 



d7 s ■ ks = 2tt cos 9 



V --\(l + y-x-5 2 )+ 
x 2 



The coefficients ctj have been determined for both decay channels. We summarize our 
results in Table || for a decaying t with positive helicity. Expressions for the other helicity 
state, or for the charge conjugate i decay, are obtained by an overall change of sign. Our 
results agree with those presented in [23], except that their expressions for the scalar case 
do not contain the factor (A 2 — B 2 )/(A 2 + B 2 )f(£, A, B). This factor contains all the 
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Table 2: Spin analyzing coefficients Oj for different decay products in SM and scalar decay of a 
polarized top quark (t| — * bW + /H + — > bl + vi) 1 or equivalently (if — > — > 6du). 

Analyzing Decay channel 

particle VK + (w = m^/m 2 ) H + (£ = m 2 H+ /mf) 

, . (l-w)(l-llo;-2a; 2 )-12u; 2 ma; 1 - £ 2 + 2£ ln£A 2 - B\ 
Vl {U) (l-^) 2 (l + 2^) (i_ 2 A 2 + m f^ B ) 



dependence of the spin analyzing power on the Lorentz structure of the coupling. We note 
further that the expression we obtain for an is in agreement with that of [44], where also 
0(a s ) corrections to this quantity are given. The inclusion of NLO corrections does not 
modify the tan/3 dependence of the ati, even if the numerical values are slightly altered. 
Since we aim to compare the analytic results to a LO Monte Carlo simulation, we use only 
the LO results for a, throughout this work. 

Since the top quark spins are not directly observable themselves, what will be acces- 
sible are quantities constructed only from the final state momenta. The most direct such 
being the doubly differential distributions of the same type as in Equation ( 3, 13] ), but now 
involving two particles one from each decaying top quark. The helicity angles 9i and 

9j are then calculated in the rest systems of the respective parents 2 . The most general 
expression of this type is 

1 d 2 N 1 , 

' 1 + Pi a, cos 9i + Piotj cos 9j + Ca^ay cos 9i cos 9j ) (3.15) 



iVdcos 9idcos9j 

where P\{P2) measures the degrees of transverse polarization of the t(t). C is the correlation 
parameter discussed in Section |2[ In leading order QCD, with the spin quantized in the 
helicity basis, P\ = P2 = by parity invariance. Instead of ( 3.15| ) the simpler distribution 

1 d 2 iV If \ 

— -r — — : — = - 1 + CdiOij COS 9; COS 9j (3.16) 

N dcos9idcos9j 4 V / K J 

is therefore expected to obtain. 

It is also possible to form one-dimensional distributions, e.g. as studied in [29]. If we 
define the angle B%j between the vectors i and j as in the previous distribution, we get 

1 dJV 1 . 

(1 + T>a%aj cos 9ij) . (3-17) 



iVdcos% 2 



2 We use the convention of performing rotation-free boosts from the ti CM system to define the orientation 
of the t(t) rest systems. Alternatively, one could perform a rotation-free boost directly from the hadronic 
CM system. 
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Figure 5: Influence of the threshold factor on the spin analyzing coefficients. The dashed curves 
show the full expression for an while the solid curves show an in the approximation /(£, A, B) = 1 



Here the coefficient T> is related to C, but in general it has a different value. To determine 
C and T> from angular distributions, the relations 

9 

C = (cos 0i cos Oj) (3.18) 

V = — (cob da) (3.19) 
cuaj 

can be used. Conversely, when the correlation coefficients are known, these relations can 
be used to determine the product c^ay. With a leading order C = 0.319, we obtain the 
corresponding T> = —0.216. 

3.5 Analysis in 2HDM (II) 

As an illustrative example of the differences between the SM and a new scalar decay, let 
us consider in some detail the results for a 2HDM (II) with couplings from Table [I]. 
In all the following, we shall fix the top mass to m t = 170.9 GeV [47]. Starting by 
analyzing the threshold region m H + ~ mt — m^, Figure |5| shows an as a function of mjj+ 
for different values of tan j3. We see that the threshold suppression becomes significant as 
m H+ approaches the kinematic limit. However, for very large, very small or intermediate 
(~ 7) tan /3 values, we infer the threshold correction to be less than 10% also for m H + = 160 
GeV. In these regions of parameter space, the threshold factor can be effectively ignored. 
Note that this argument is not specific to an, but applies to all cej, since f(£,A,B) is a 
universal factor. 

Summarizing the results presented in Table || for the 2HDM (II), Figure |6| shows a 
numerical evaluation of the analytic expressions for all a%. The results are presented for 
two values of tan/3: one large value (tan/3 = 50), for which the efficiency in analyzing the 
spin is optimum, and one intermediate value tan/3 = \J mt /m^ ~ 8 (for mt, = 3.2 GeV 
with mjj+ = 100 GeV), where all the sensitivity to analyze the top spin vanishes. This 
value corresponds to a purely scalar coupling, thus to an isotropic decay. For tan (3 < 8, 
all cc, acquire an extra minus sign compared to tan (3 > 8. This corresponds to a shift 
from predominantly right-chiral to left-chiral coupling. We see also, that the efficiency to 
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Figure 6: Spin analyzing coefficients cc, for top decay within SM (solid black), and in the 2HDM 
(II) (dashed red), presented for tan/3 = 50 (left) and tan/3 = 8 (right). For low tan/3 ~ 1, the 
Higgs results presented in the left plot are valid with an overall change of sign. 
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Figure 7: The solid red curve shows at for the 2HDM (II) , whereas the dashed curves correspond 
to the SUSY-corrected model (II) with tb = —e' t = ±0.01. The mass m H + = 100 GeV was used. 
For comparison, the black line shows ab for the SM decay. 

analyze the top spin is not highest using the charged lepton, as is shown to be the case 
for the SM. Instead the most efficient probe is either the Higgs momentum itself or the 
associated b quark. This is easily understood; since the itself does not carry any spin, 
the top spin information can only be transferred to the angular distributions of the b/H + . 
In a vector decay, parts of this information go into the different polarization states of the 
W, as discussed above. 

Figure [7| displays a& as a function of tan (3 for a fixed m#+ = 100 GeV. It illustrates 
clearly the transition from = — 1 analyzing power at small tan/3, to a\, = 1 associated 
with a right handed coupling for large tan (5. Note also, that for tan [3 = */ \f2"^- ~ 6, the b 
quark coefficient a& of the Higgs events mimics that of the SM decay. The only dependence 
in this relation on mjj+ enters through the running mass m^. No corresponding value exists 
where the a; are equal, since the Higgs value is bounded by kinematics to \ol\ \ < 0.5, whereas 
the SM value is always a\ = 1. 

To compare with the case when tan /3-enhanced SUSY corrections to the charged Higgs 
couplings are included, we give in Figure ^ also the values of for the 2HDM (II). Rather 
than to calculate the corrections for a specific SUSY model spectrum, we parametrize them 
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Figure 8: Upper row shows doubly differential distributions in cosf^j (denned in the text) of the 
decay products d and r from ti decaying through W ± W T (left) and H ± W T (center, right). Results 
obtained at ME level for m#+ = 80 GeV, tan/3 = 50 (center) and tan/3 = 1 (right). The bottom 
panel lines correspond to projected angular distributions with cos 0d > (solid) and cos 6^ < 
(dashed). 



in terms of the parameters | et, | < 0.01 and \e' t \ < 0.01. These are reasonable maximum 
values [55], which correspond roughly to a s /37r as discussed above. Figure clearly shows 
that even though these corrections are enhanced by tan/3 in the couplings, they have only 
a small effect on the ratio that enters the spin analyzing coefficients. In fact, the largest 
correction is obtained not in the high tan (3 limit, but in the transition region around 
tan (3 = 8 - 20. Given the observed smallness of the SUSY effects on the spin analyzing 
coefficients, it is acceptable to apply the results from the 2HDM (II) without modification, 
both in the high and in the low tan/3 regimes. However, total rates for t — > bH + are of 
course affected by the differences. 

A few words are to be said also about the 0(a s ) corrections calculated in [44]. Inclusion 
of these effects leads to modifications of the spin analyzing coefficients in a fashion very 
similar to the tan/3 enhanced SUSY corrections discussed above. These corrections are 
also largest in the intermediate tan/3 region, where they can reach 20% in magnitude. In 
the large and small tan /3 limits, the NLO corrections have negligible impact. Since our 
results are most interesting in these limits, we will show plots for tan /3 = 50 and tan (3 = 1. 
Interpolation to the intermediate range should then be performed with care, remembering 
the higher order corrections. 

As a final result for the 2HDM (II) on matrix element level, we show the differen- 
tial distributions in (cos 9i, cos Oj) and cos described by Equations (3.16) and ( 3.17] ) 
respectively. The top row in Figure || shows lepton-lepton correlations, where the lepton 
(d quark) from a W decay is correlated with another lepton (in this case a r) from 
or decay from the opposite side of the event. In the absence of spin correlations, this 
distribution is expected to be flat. From left to right, Figure gives the results for the 
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Figure 9: Angular distributions denned similarly to those in Figure |[ with the r replaced by the 
associated b quark originating from the same t(t) decay. Again m#+ — 80 GeV and the ordering is 
SM (left), tan/3 = 50 (center), and tan/3 = 1 (right). 



SM, for the 2HDM (II) with tan/3 = 50, and similarly with tan/3 = 1. In the bottom row, 
distributions in one of the angles are given. These are obtained using the other angle to 
determine the parent spin by applying the projections cos 6^ > and cos^ < to the 2D 
distributions. Figure ^ illustrates the points we have previously made about the 2HDM 
(II), namely that i) the lepton-lepton correlation is more efficient in the SM since a.\ ^ 1 
in the Higgs case, ii) there is a change of sign in going from high to low values of tan/3, 
and hi) the low tan (3 case is more SM-like. 

Going to Figure ^, the r lepton has been replaced with the b quark associated with 
the same (Higgs) side of the event. The three distributions are defined similarly to those 
in Figure ||, as are the projections. From Figure ||, the increased spin analyzing efficiency 
when using the associated b quark in the Higgs case is evident. 



The last figure to discuss here is Figure 1C, which shows distributions in cos 9™ between 
different particles in the final state. An uncorrelated sample corresponds to a flat line 
distribution. The specific combinations chosen are (r,d), (b,d), (b,b), and (t, b), where 
the first particle originates from H^/W^ decay, while the second always comes from the 
opposite side VF ± . The distributions (H + , b) and (H + , d) are equivalent to the distributions 
(b, b) and (6, d) displayed in Figure [To] with a change of sign in the spin analyzing coefficient 
(au = —db in Eq. J3.17 )). Hence they would not contribute any additional information. 

By combining in the same plots for the SM with results from the 2HDM, we see directly 
which correlations are more efficient in the two cases. In agreement with previous results, 
this is again the rd distribution for the SM, and the bb distribution for the 2HDM (II), 
illustrating the universal dependence of the spin correlation effects on c^ay. We have also 
found, that when cuts are applied at the parton level, the distributions shown in Figure 10 
are not affected nearly as much as those presented in Figures ^ and |9[ Using the cos 9{j 
variables to study ti spin correlations is therefore advantageous to avoid the problems with 
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Figure 10: Distributions at matrix element level of cos Oij for different final state particles 
The dotted line is the expected SM result without tt spin correlations. The SM result is shown 
in black, while the 2HDM (11) is presented for tan/3 = 1 (long-dashed blue), and tan/3 = 50 
(short-dashed red). All results are given for m H + = 80 GeV. 



cuts discussed in [23]. We will return to this discussion below. 



4. Monte Carlo Simulations 

Up to this point, the results we have presented were obtained directly from matrix elements. 
To give a more realistic assessment of the prospects to observe any of these spin effects in a 
collider experiment, complete hadron- level events must be considered. We do this using a 
Monte Carlo (MC) approach. The framework is again the 2HDM (II) because of its special 
status as the minimal Higgs model compatible with supersymmetry. 

The production of ti is treated completely within the SM. On the decay end, we study 
in parallel the situations when either both decays occur within the SM, or when one of 
the two top quarks decays through W and the other through H^ 1 . 3 Being an interesting 
process in its own right, the SM decay of ti is also the main irreducible background in 
searches for light [14]. Over nearly the full range of tan/3, if ± decays preferentially to 
the heaviest lepton available, that is H + — > t + v t . We therefore restrict ourselves to this 
decay channel. Consequently, for the SM events, we demand one of the two W bosons to 
decay through this mode which has BTZ(W ± — ► t^u t ) = 0.1125. The r lepton subsequently 

3 Of course, there is also the possibility of having both tops in an event decay through . However, since 
this mode gives a different final state, and since the effect is sub-leading, we will neglect this contribution 
here. 
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decays hadronically producing a r jet. For the other W, which is always present in the 
event, we consider the hadronic decay to allow for hadronic top reconstruction. We show 
SM and 2HDM results separately, keeping in mind that BTZ(t — * bH + ,H + — > t + u t ) and 
BlZ{t —* 6V1 /+ ,V1 /+ — * t + v t ) are of similar magnitude for the tan/3 and mjj+ regions of 
interest. 

To incorporate helicity information throughout the whole process, it is necessary to 
use a MC generator which can treat the full 2^6 matrix element. This is provided by 



MadGraph/MadEvent 4.1.10 [63] for which we implemented the model (3.2). Using this 
program has the additional advantage that the ME is generated without using the narrow 
width approximation. To verify the treatment of spin information, control samples were 
generated using the specialized ti generator TopReX 4.11 [64]. For the observables we 
analyze, results from both programs were found to be in good agreement (after correcting 
the partial width T(t — ► bH + ) in TopReX). 

From the 2 — ► 6 matrix element generated by MadEvent, full events are obtained by 
applying Pythia 6.409 [65] for parton showering and hadronization. To treat properly the 
spin information in tau decays, which was previously demonstrated to be important in 
searches [66,67], Tauola [68] is invoked. The underlying event is modeled using the Pythia 
default "old" model based on multiple parton-parton interactions. The default parame- 
ters of this model are tuned to Tevatron minimum-bias data [69] and provide reasonable 
estimates for extrapolation to the LHC energy. 

4.1 Event Reconstruction 

For jet reconstruction we use the Fast Jet implementation [70] of the longitudinally invariant 
kj_ algorithm [71-73] for hadron colliders. The jet clustering uses i?-scheme recombination 
based on the distance measure dij = mm{k 2 ±i , kl^ARfj/R 2 , with AR = ^(Arj) 2 + (Ac/>) 2 
and R 2 = 1. We use the algorithm in exclusive mode, meaning that particles are clustered 
with the beam when the beam-particle distance diB is smaller than the distance dij to 
any jet candidate. Not all particles will therefore end up in a jet. Furthermore, we only 
take particles with \rj\ < 5 into account to reflect the detector acceptance region. For the 
minimum jet separation measure, above which no further clustering takes place, the value 
d C ut = 400 (GeV) 2 is used. Choosing this value gives a jet multiplicity which peaks at the 
value expected from the matrix element. 

We also implement a simplistic notion of flavor tagging where jets are tagged as b jets 
or t jets by comparing to MC truth information. A candidate jet is tagged whenever the 
distance AR to a true b quark (or a r) is less than 0.4. In addition, for the jet to be tagged, 
it is required that it has \rj\ < 2.5. Apart from these criteria, no further efficiency factor is 
used in the flavor tagging. 

Events are selected for analysis based on their overall topology. The characteristic 
signature, which must be fulfilled by our signal events following jet reconstruction, is the 
presence of exactly two b jets, one r jet and at least two additional, untagged, jets. For this 
type of potential ti events, the hadronically decaying W is reconstructed by combining two 
light jets, minimizing the mass-square difference Am^, = m 2 - — m^. For W candidates 
within a certain mass range Amw < 10 GeV of the true W mass, further recombination 
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with one of two b jets is performed to reconstruct a t(i) candidate. On this candidate, a 
similar cut as for the W on Amf = fn 2 -- h — mf is applied to asses the overall goodness of 
the reconstruction. We keep events which have Amj < 15 GeV. For the surviving events, 
one side has then been fully reconstructed. The remaining b jet, not used in the top 
reconstruction, can therefore be associated with the t — ► bH^/W^ — > bT^v T decay. 

To show the statistical efficiency, we give in Table ||] the total number of MC events 
generated for each model. Since the branching ratios to reach the requested final state are 
largely different for the SM and the 2HDM, this nearly constant number of events actually 
corresponds to quite different integrated luminosities, as also given in the table. Note 
that, for these low Higgs masses, the luminosity corresponding to this number of events 
is highest for the SM sample. Table |3] also lists the total number of events left in each 
case after reconstruction. These numbers indicate the level of statistical uncertainty in the 
jet-level distributions we show below. 

Due to the presence of at least two neutrinos in the final state from the decay of the 
(or W^) into t^v t , the longitudinal momentum of the tit) on that side of the event 
cannot be reconstructed. Thus its rest system is not accessible, neither is the overall CM 
frame. This fact effectively prevents the direct experimental use of distributions such as 
(3.16) or (3.17) to establish the presence of a t —* bH + channel. This is in contrast to the 
SM case where the semi-leptonic (e,fj) channels could be very useful in establishing the 
existence of ti spin correlations. In the following sections, we will focus on how to address 
this somewhat discouraging situation. 



4.2 Distributions in A<ft 

We will first investigate if there are interesting observables defined directly in the laboratory 
frame. An early study [19] established that distributions in azimuthal angle A<ft(l + l~), 
with the l + l~ originating from two leptonic decays of a ti pair, are sensitive to ti spin 
correlations. These results were recently confirmed in a NLO shower MC study [30], where 
also other spin-sensitive observables defined directly from lab quantities were considered. 
Their conclusion was that, for all the variables included in the study (see [30] for the 
complete list), the spin effects are only visible in the distribution A(p(l + l~). 



Table 3: Number of events generated at MC level, and events left after reconstruction. Integrated 
luminosities, corresponding to the generated event samples, arc calculated using the NLO ti cross 
section. The luminosity varies due to the different branching ratios for the final state rv T jj. 
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Figure 11: Azimuthal angle between b jet from W decay and opposite side b- or r jet. The 
results are shown for the SM ti (black), or t — ► bH + decay with tan/3 = 50 (short-dashed red) and 
tan/3 = 1 (long-dashed blue). For comparison, the distribution obtained in the SM without spin 
correlations is shown as a dotted black line. The mass m H + = 80 GeV was used. 



Since we require hadronic reconstruction of a W boson and a top quark, we are effec- 
tively limited to studying the two distributions A(j)(bb) and A<^(r6). From what we know 
about the spin analyzing efficiencies, we could expect that, at least for the Higgs case, 
these variables could show some effects which did not show in the SM study. The resulting 



azimuthal distributions are shown in Figure 11. In addition to showing the SM and the 



2HDM (II) for tan/3 = 50, and tan/3 = 1, we show also this time the SM without any spin 
correlations to illustrate that this distribution is not a flat line. From the figure, in the bb 
distribution, it can be seen again that the low tan (3 case is more SM-like and that the high 
tan/3 case gives the most significant deviation from the SM distribution. In the rb distri- 
bution, the curves are hardly separable. It therefore remains to be seen if these variables 
could be used in any way to discriminate between the SM and new physics contributions. 

Finally we note that these distributions are quite sensitive to changes in the Higgs 
mass. Larger values of A</> are favored with increasing m H +. The kinematic effects must 
therefore be under good control if the A(f> distributions should be used to determine tan (3. 

4.3 Brandenburg observables 

We now turn our attention from the laboratory frame observables to what can be done 
with half the event fully reconstructed. This has been discussed by Brandenburg who 
introduced [26] the observables 

Ox = • k 2 

2 = (k;-z)(k 2 -2) (4.1) 

3 = (k;-p)(k a -p), 

where an asterisk denotes momenta in the rest frame of the parent top quark. The beam 
direction is given by z, and k 2 is the laboratory momentum of the lepton from the top 
decay not associated with ki. Each of the observables 0\ - O3 are cleverly constructed 
to have mean value strictly equal to zero in the absence of ti spin correlations. We first 
evaluated 0\ - O3 at parton-level, similarly to what is done in [26], but also including 
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the decay through H^ 1 . The resulting mean values obtained from the distributions in 
these observables are presented in Table ^. Their standard deviations are AO\ = 0.58, 
A02 = 0.42 and AO3 = 0.39 for all three cases. Comparing our SM values with those 
given in the original study, we find that there is an overall agreement. For the Higgs case, 
we observe the same sign differences between the low and high tan j3 regimes as previously 
seen for the correlation plots. We also note that the Higgs values are generally smaller 
than their SM counterparts. 

Unfortunately, when evaluated at the jet level, including the cuts described above, we 
find that the expectation values for all three Brandenburg observables are consistent with 
zero. 



4.4 Transverse correlations 



As an alternative to the full reconstruction of the rest frames for t and t, and the half recon- 
struction performed when calculating the Brandenburg observables, we consider transverse 
reconstruction. This is possible in the hadronic channel where one top quark is recon- 
structed using the methods discussed above, while the transverse momentum of the other 
top can be obtained from the sum p±j = p± t b + p_\_ :T + p_i_, m iss- The sum involves the b jet 
not used for top reconstruction in the first step. 

Following the reconstruction, we perform boosts of particles into the transverse 
rest frames of t(t), and consider distributions in this frame of azimuthal angles (Acfti, A(f)j) 
to the p± axes of t(t). It is then possible to form both doubly differential distributions 
similar to ( [3.16[) , or one-dimensional distributions similar to (|3. 17 ) . To obtain a distribution 
which like ( 3.17| ) is robust with respect to cuts, we have looked specifically at distributions 
in Acpi — Acf)j. In analogy with the construction of CM variables, we expect a distribution 
of the form 

- - 1 + V'ctiCtj cos(A0j — Aifij] 



Ndcos(A4>i - A(f>j 



(4.2) 



This relation involves yet another spin correlation coefficient T>'. To illustrate how similar 



this new observable is to cos 9{j, which was previously shown in Figure 10, Figure 12 shows 
the matrix element results on cos(A(/>j — Acj)j) for the same particle combinations. The 
results are presented for the SM, and for the charged Higgs decay with tan/3 = 50 and 
tan (5 = 1. From comparing the two figures, we note that T>' ~ 0.9P, signaling only a slight 
loss of correlation going to the transverse projection. 



Table 4: Mean value of the Brandenburg observables in percent. Evaluated for SM (W^VJ^T) and 
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Figure 12: Distributions at matrix element level of cos(A</>j — A</>j) for different final state particles 
Dotted line show results in the SM without spin correlations. SM results with correlations 
are shown in solid black. The 2HDM (II) is shown for tan/? = I (dashed blue), and tan/3 = 50 
(dotted red). All results were obtained for m H + = 80 GeV. 



These transverse correlations also maintain their properties when going to the level 



of jets, as demonstrated in Figure 13. That only two diagrams are shown in Figure 13 



and not even the ones from Figure 12 with best separation between the different cases, is 
because the efficiency in experimentally separating a d-type jet from a -u-type jet is only 
61% [23]. Since this separation step is required in order to use the jets from the hadronic 
W decay in the analysis, we exclude the (r, d) and (b, d) combinations, which both require 
this information. We remind the reader that the distribution (H + ,b) contains the same 
information as that for (6, b). This is true also at jet level, since the momenta of H + and b 
are still anti-parallel in the reconstructed top rest frame. Figure 13 shows good separation 
between the SM and the 2HDM (II) for tan/? = 50, while the small tan/3 case is more 
difficult to separate from the SM. Since this small separation is partly intrinsic, as can be 
seen from the lower panels in Figure 12, these variables shows most promise in the high 
tan/3 regime. 



We have studied the dependence of the distributions in Figure 13 on m^+. From 



their a dependence, the coupling sensitivity in the rb correlation is expected to degrade 
for higher mjj+, while it should not change for bb. This is indeed what is observed. We 
also find additional distortion in both distributions when mjj+ > 130 GeV. For bb, values 
cos(A<^, — A^) < are favored, while for the rb distribution instead cos(A^ T — A0g) > 
shows a surplus of events. However, due to the small values of a T for , the rb distribution 
is hardly interesting for such masses. This additional angular dependence, which is not 
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Figure 13: Distributions at jet level of cos(A0i — A<j>j) for different final state particles 
The different contours correspond to the SM with (solid black) , and without (black dotted) , tt spin 
correlations. There is also the 2HDM (II) for tan/3 = I (dashed blue), and for tan/3 = 50 (dotted 
red). All results were obtained for = 80 GeV. 



observed at parton level, we attribute to the b jet reconstruction. For higher m H +, the b 
quarks from t — > bH + becomes softer, and eventually falls below the jet measure ri cu t in 
p±. These events should normally not pass for further analysis, if not the underlying event 
provided the necessary energy. When these jets are dominated by the underlying event, 
their directions are smeared, causing the t rest frame to be poorly reconstructed. The 
result is an additional distortion in the angular distributions compared to the ME results. 



5. Summary and conclusions 

Top quarks produced in pairs at the LHC should experience strong correlations among 
their spin projections in the helicity basis. This entanglement makes it a favorable system 
for studying the Lorentz structure of the couplings involved in top quark production and 
decay. We have considered here the effects on spin dependent angular observables of new 
physics in top quark decays. Specifically, we have discussed top decay through a charged 
Higgs boson. 

Within a general CP-conserving 2HDM, we obtain results on the spin analyzing coef- 
ficients. They determine the sensitivity of a given decay product to the spin projection of 
the decaying particle. Unlike the SM decay t — ► bW + — > bl + vi, where the charged lepton 
(or d quark jet) is the most powerful spin analyzer, the scalar case has the associated b 
quark, or the momentum itself, as the most efficient spin analyzers. As for the lepton 
in the SM case, this efficiency can reach unity. In the 2HDM (II), this is found to be the 
case for very large (and very small) values of tan (3. At the level of matrix elements, the 
modification to the distributions of certain angular observables, induced by the presence of 
an H^, could therefore be similar in magnitude to the SM spin effects themselves. Whether 
this is actually the case depends on the parameters of the 2HDM. Since all spin analyzing 
coefficients associated with t — ► bH^ — * br^v T depend universally on the Lorentz structure 
of the tbH + vertex, an angular analysis of the decay products could provide an additional 
handle on the couplings. This includes the possibility of determining the effective 
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value of tan (3 in the MSSM, since we find that tan (3 enhanced SUSY corrections to the 
spin correlation observables are small. 

The preference of the to decay into t^v t prevents reconstruction of longitudinal 
momenta. Thus the t(t) rest frame spin analysis, which can be performed in the SM, 
cannot be directly applied to the new physics case. Encouraged by the significant spin 
effects observed at matrix element level, we have therefore studied several longitudinally 
boost-invariant observables. 

We find that the distribution in azimuthal angle between the b jet associated with 
t — ► bH + , and the b jet from the opposite side t — ► bW decay, shows some sensitivity to the 
Lorentz structure of the coupling. Although the effects on this observable are small, 
they are directly measurable and therefore shows promise for further experimental study. 

Even more promising are (A^, A(f>j) correlations in the transverse rest frames of tit). 
This observable is also directly measurable, and it is found to be robust with respect to 
phase-space cuts. The highest sensitivity is also here obtained for the bb correlation. This 
distribution should be particularly useful in the highly interesting case with large tan (3, 
as it lifts the degeneracy between large and small tan/3 values present in the inclusive 
B1Z(t — > bH + ) measurement. It should be investigated if this observable could be useful 
also for studying ti spin effects within the SM. In that case, it should be possible to use 
lb distributions where I = e,[i. Another question is if the correlations in bb or lb obtained 
from the semi-leptonic ti decay in the SM could be used as control samples. 

As illustrated by our MC samples, statistics is not the limiting factor in an analysis of 
this type due to the large cross section for production at the LHC. We therefore conclude 
that further investigation of spin effects in top physics beyond the Standard Model should 
be worthwhile. 
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A. Massive Spinor formalism 

In this appendix, we discuss the formalism used to describe polarized massive fermions. 
Throughout the discussion, we try to adhere to the conventions of [74] which contains an 
excellent introduction to spin formalism. 

For massive spin- 1/2 particles, the projection relations 




(A.l) 
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obtained by summing over the possible helicity states A are very well-known, as are the 
more general projector relations 



2 

with a time-like spin four- vector defined as 



u(p, X)u(p, A) = i (1 + 75^) {f + m) 
v(p, \)v(p, A) = - (1 + 75^) (f - m) 



(A.2) 



2A(M^p). (A.3) 



m m 



In the rest frame of the massive particle = 2A(0,p), whereas in the high-energy limit 
= 2\p fJ "/m. As a further generalization of the projector technique, Bouchiat and Michel 
[75] introduced a projection relation for the more general product of two Dirac spinors of 
mass m with arbitrary helicities. For this we need to introduce three spin vectors s° that 
are mutually orthonormal, and in addition orthogonal to p/m in the sense: 



p ■ s a = 

(A.4) 



s a -s b = -5 ab 



s aa__ PjJPu 



The most direct choice for an explicit spin basis of this type is perhaps, in a coordinate 
system where p is in the positive z-direction, to use 



(A.5) 





(0,x) 




(o,y) 


s 3 M = 


V m m 



With this choice, s 3>M differs from the positive helicity spin vector ( A.3 ) only by a factor 
2A. The helicity spinors can be shown to satisfy 

j 5 fu(p,\') = al x ,u(p,X) 
T J a v(p,\') = a a yx v(p,\) 

where the rhs is summed over A, and the a a are Pauli matrices where the first (second) row 
and column correspond to A = 1/2 (A = —1/2). Thus cr 3 = 2 XI. From these identities, the 
Bouchiat-Michel relations 

u(p, X')u(p, A) = - (5 xy + 7 5 f a a xy ) {f + m) 

I (A.7) 
v( P , X')v(p, A) = - (<5aa' + 7 5 ^Vva) [f-m). 

follow. Letting A' = A (no sum over A), and noting that s\ x ^ only for i = 3 with the 
spin basis ( |A.5| ), the relations ( |A.2|) are obtained as a special case of the Bouchiat-Michel 
formulae. 
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